Nonuniform mixed-parity superfluid state in Fermi gases 
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We study the effects of dipole interaction on the superfluidity in a homogeneous Fermi gas with 
population imbalance. We show that the Larkin-Ovchinnikov-Fulde-Ferrell phase is replaced by 
another nonuniform superfluid phase, in which the order parameter has a nonzero triplet component 
induced by the dipole interaction. 
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Bardeen-Cooper-Schrieffer (BCS) theory is the stan- 
dard model of superconductivity and superfluidity in 
fermionic systems. The central concept of the BCS model 
is that of a Cooper pair composed of fermions with oppo- 
site momenta and spins. If there is a mismatch between 
the Fermi surfaces of spin- up and spin-down fermions, 
e.g. because of the Zeeman splitting in magnetic field, 
then the formation of the Cooper pairs costs energy, 
which results in the suppression of the critical temper- 
ature. It was shown by Larkin and Ovchinnikov and 
Fulde and Ferrell 0] (LOFF) that the pair-breaking ef- 
fect of the Fermi surface splitting can be reduced if the 
Cooper pairs have a nonzero center-of-mass momentum. 
The resulting nonuniform superconducting state turns 
out to be more favorable at low temperatures than the 
uniform BCS state. 

Due to the sensitivity of the nonuniform state to dis- 
order and the orbital effects, the attempts to find it in 
superconductors have been unsuccessful (at least until 
the recent observation of the features consistent with the 
LOFF state in the phase diagram of the heavy-fermion 
compound CeCoIn 5 The recent surge of interest to 
the nonuniform superconducting and superfluid states, 
see e.g. Refs. 0, 0, 0> S H> @> nas been stimulated 
by the experimental progress in ultracold atomic Fermi 
gases, such as 6 Li and 40 K. The Fermi-surface splitting 
in these systems is due to the populations of atoms in 
different hyperfine states being unequal 0, . When 
a pairing interaction between the fermions is turned on, 
the system becomes formally equivalent to a neutral and 
perfectly clean superconductor in a Zeeman field. 

In this Letter we consider the effects of the long-range 
and anisotropic interaction between atomic dipole mo- 
ments in a Fermi gas with population imbalance. The 
magnetic dipole interaction is quite small for alkali atoms 
but can be considerably enhanced to become experimen- 
tally observable for atoms with large dipole moments ^3] • 
Another possibility is to control the magnitude of elec- 
tric dipole moments by applying external electric field, as 
suggested in Ref. 01 ■ I n either case, we assume that the 
dipole interaction represents a correction to the dominant 
isotropic s-wave pairing. We will show that increasing 
the population imbalance produces a novel nonuniform 



mixed-parity (NMP) superfluid state, which has higher 
critical temperature than the LOFF state. 

We consider a homogeneous Fermi gas consisting of two 
different species of atoms of equal mass m. For instance, 
in 6 Li one can have a mixture of the hyperfine states 
M = \m„ = 1/2, mi = 1) and |-) = \m s = 1/2, m 4 = 0) 
|14j . where m s and rrii are the electron and nuclear spin 
projections respectively. The free-particle Hamiltonian 
has the form H = J2k(^S a p - ha s , a ff)cl a Ck^, where 
£k = (k 2 — hp) /2m, Uf is the Fermi momentum, h is 
the Fermi surface splitting due to unequal particle con- 
centrations in the states labelled by a, = ±, and 03 
is the Pauli matrix (we set h = ks = 1). We assume 
that the atomic dipoles are fully polarized, e.g. by an 
external magnetic field along the z-axis: fi = (A4Z. The 
interaction between two atoms at distance R is 



U(TL) = -A <5(R) + A 



1 - 3Rj 
R 3 



(1) 



where the first term describes a short-range attraction 
with the coupling constant Ao = 47r|a|/m > 0, and the 
second term is the dipole interaction. The s-wave scat- 
tering length a < depends on the field, diverging in 
the strong-coupling regime near the Feshbach resonance. 
The effects of the dipole interaction JIJ have been exten- 
sively studied for Bose gases, see e.g. Ref. and the 
references therein. 

In the BCS regime the pairing Hamiltonian, which 
takes into account the possibility of the Cooper pairs 
having a nonzero momentum q, can be written in the 
following form: 
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(2) 



where V(k,k') = -A + 47r/4(*k_r ^) 2 /|k - k'| 2 is the 
Fourier transform of Eq. JIJ [lfj|. The matrix T has 

the following nonzero elements: = |_ = 

r = r = 1, which is due to the fact that 

the pairing interaction depends only on the total particle 
densities n(r) = n + (r) + ri_(r). 
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Treating the Hamiltonian J5J in the mean-field ap- 
proximation, one arrives at the self-consistency equation 
which contains a divergent momentum integral due to 
the absence of ultraviolet cutoff in V(k, k'). This di- 
vergence is removed by the renormalization of the in- 
teraction potential |l7|. in which the bare interaction is 
replaced by the two-particle scattering amplitude, and 
an effective cutoff cu c , of the order of the Fermi energy 
e F , is introduced. We use a BCS-like sharp cutoff and 
represent the pairing interaction in a factorized form: 

v(k,k') = -e(\&\ -w e )e(\&\ -u^E.A^k^k'). 

The symmetry factors 4>i(k) are the eigenfunctions and 
the coupling constants — A^ the eigenvalues of the opera- 
tor V defined by the kernel V(kpk, kpk.'). The symmetry 
factors are even (odd) for the singlet (triplet) pairing and 
satisfy the orthonormality condition: (0*(k)0j(k)) =<5,j, 
where the angular brackets denote the average over the 
spherical Fermi surface. 

The effect of the dipole interaction on superfluidity is 
two-fold: in addition to breaking the rotational symme- 
try of the system and introducing gap anisotropy, it also 
leads to the possibility of triplet pairing. We include 
one singlet and one triplet pairing channel with the high- 
est critical temperatures. One can show that the lowest 
eigenvalues of V in both cases are achieved for the eigen- 
functions that do not depend on the azimuthal angle: 
0(k) = f(cos8), where the function /(s) satisfies the in- 
tegral equation J^ 1 ds'K(s,s')f(s') — — A/(s), with the 
kernel K(s, s') — — A /2 + n^\s — s'\. In the singlet case 
the solution is f(s) cx cos(ks), while in the triplet case 
f(s) tx sin(/is), where k — y/2irn^/A. For the singlet 
pairing, the normalized symmetry factor has the form 



<Mk) 



■ cos(/t s cos ( 



where k s satisfies the equation 



2^ 



A - 2-rrfi 2 , ' 



(3) 



(4) 



from which one obtains X s = 2tt^i 2 ,/k 2 . In the limit 
of weak dipole interaction, /j 2 , <C Ao, we find </> s (k) ~ 
1 + n 2 {l — 3s 2 )/6, and A s ~ A — 27r^, i.e. the dipole 
interaction leads to a downward renormalization of the 
coupling constant in the singlet channel. For the triplet 
pairing, one can show that the lowest eigenvalue corre- 
sponds to k — 7r/2, so that 



t(k) = \/2sm(^cos6») 



The singlet and triplet contributions to the pairing 
Hamiltonian J5J can now be explicitly separated: 



Hi, 



-A a £Bt(q)B fl (q) - At^B^IMq), (6) 



where Bj(q) = (1/2) £ k ^(k)(iCT 2 ) a/3 4 +qAa u_ k+qA/3 

and Bj(q) = (1/2) £ k 0t(k)(io-<T 2 ) Q/ 3cj c + q /2,a c -k+q/2, 1 8 
are the pair creation operators. The representation © 
emphasizes the full invariance of the pairing interac- 
tion with respect to rotations in the "pseudospin" space 
spanned by the states |+) and |— ). Decoupling Hi nt 
yields the gap function 

A Q/3 (k,r) = (icr 2 )Q/3^(r)0 s (k)-(-(icrcr 2 ) ct/3 d(r)0 t (k), (7) 

where the complex scalar if) and the complex vector d 
are the order parameters of the singlet and the triplet 
pairing respectively. 

In the absence of population imbalance h — 0, and the 
critical temperatures for the singlet and triplet pairing 
are given by the standard BCS expressions: T a = = 
(2e c /7r)w c e -1 / A ' FAa , where a = s,t and Np is the density 
of states at the Fermi surface per one hyperfine state (all 
three components of d have the same critical tempera- 
ture). Of the most interest to us is the limit of weak 
dipole interaction, in which T s 3> T t . 

The phase diagram at h can be obtained from 
the free energy, which is represented as an expansion 
in powers of the order parameter components: T = 
J~o + JF^^d] + ■•• {Fq is the free energy in the nor- 
mal state). To find the critical temperature T c (h) of the 
second-order phase transition into the superfluid state, 
or inversely the critical Fermi-surface splitting h c (T), it 
is sufficient to keep only the quadratic terms in T . Using 
Eq. © we find that the contributions to Ti from ip and 
d Zl which describe the pairing between fermions of dif- 
ferent species, are decoupled from d± = (d x ± id y )/^/2, 
which correspond to the intra-species pairing: 
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(8) 



The coefficients in this expression have the following 
form: A ab = A ba = N F [S ab In (T/T a ) + I ab ], A± = 
N F []xi(T/T t ) + I a (h = 0)±SI], 



(5) lab = ( <Mk)<Mk)Rc * 



1 
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and A t = (8/tt)^ 2 ,. Note that the coupling constant and 
the symmetry factor in the triplet channel are not sensi- 
tive to the value of An and therefore are the same as in 
the pure dipolar case |18| . 



where Sit (k) are defined by Eqs. © and (JSJl, $>(x) is 
the digamma function, W = i>Fkq — 2h, vf is the Fermi 
velocity, and SI cx N' F h is proportional to the difference 
between the densities of states at the Fermi levels for 
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the two fermionic species. Since at h ^ the symmetry 
in the "pseudospin" space is reduced to rotations about 
the z-axis, different components of d appear at different 
temperatures. 

The critical temperatures for d + and d— are found from 
the equations A + = and A_ = respectively. If one ne- 
glects the band asymmetry N' F then the phase transition 
is not affected by the population imbalance, otherwise 
T c ,±(h) — T t =F O(SI), so that a nonunitary triplet state 
with d+ = 0, d- 7^ is realized. 

The effect of the Fermi-surface splitting on the other 
two components of the order parameter, ip and d z , is more 
interesting. Because of the presence of the off-diagonal 
matrix elements in A, see Eq. (0, the singlet and triplet 
pairing channels can be mixed at q ^ and ft / 0, 
producing the NMP state, in which both ip and d z are 
nonzero. The critical temperature is obtained from the 
equation det A(q) = 0, after maximization with respect 
to q. The calculation is facilitated by the observation 
that the functions Q can be expressed in terms of two 
dimensionless variables Q = VFl/2h and z — h/2nT. At 
any given < z < oo, we find 

T c {h) =T S max e I( - Q - z \ (10) 
Q 

where 

% = -\{Iss + ht +r) + I v /(/ Sfl -/ tt _ r )2 +4 /2. 

The critical band splitting is given by h c (T) = 2irzT c . If 
the maximum of X is achieved at Q = Q c , then the wave 
vector of the superfluid instability is q c = 21i c Q c /vf- 
The parameter r = \n(T s /T t ) > characterizes the rela- 
tive strength of pairing in the singlet and triplet channels. 

In the absence of dipole interaction, r — oo, the pair- 
ing is isotropic, and J(Q, z) = #(1/2) - (Re^(l/2- iz + 
izkQ)). At z < zloff — 0.30, this function has a max- 
imum at Q = 0, therefore the phase transition occurs 
into the uniform superfluid state. At z > zloff, i- e - a t 
X < Xloff — 0.56X,,, the maximum of X is achieved at 
|Q| 7^ 0, which corresponds to the LOFF state. 

Since the dipole interaction lifts the spherical degen- 
eracy of the critical temperature, the cases q 1 z and 
q || z should be considered separately. In our numerical 
analysis we use the following values of the parameters, 
appropriate for the weak dipole interaction in the BCS 
limit: (A^Ao) -1 = 1.0 and p 2 d /\o = 0.1, which gives 
T t ~ 0.02T S and r ~ 3.93. 

For q = qz, we find that at z < znmp — 0.23 the 
maximum of X is at Q = 0, therefore I st = and the 
phase transition occurs into the uniform singlet state. In 
contrast, at z > znmp, i-e. at T < Xnmp — 0.68X S , the 
critical temperature has two degenerate maxima at Q = 
±Q c z ^ 0, which corresponds to the phase transition 
into the state 

U HO) 



Ts |— r 
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FIG. 1: Critical temperature vs Fermi-surface splitting in a 
homogeneous Fermi gas with dipole interaction. The solid line 
corresponds to the transition into the NMP state with q || z, 
the dashed line corresponds to the LOFF state with q _L z, 
and the short-dashed line is the instability line of the uniform 
BCS singlet state. At the lowest temperatures, T < Tt, the 
nonunitary triplet state (TSC) is realized. 

Here p = ~I st /[ln(T c /T t ) + I u ], and the weights 771,2 of 
the two plane-wave components are determined by min- 
imizing the full nonlinear free energy J 7 , see below. For 
q _L z, there is no singlet-triplet mixing terms, and at 
T < Xloff one obtains the LOFF state with tp modu- 
lated in the equatorial plane and d z = 0. The transi- 
tion line for this state is determined by the maximum 
of X(Q,z) = —I ss , which yields the critical temperature 
slightly higher than in the isotropic LOFF case but still 
lower than in the NMP case. 

Thus we come to the conclusion that it is the NMP 
state l|ll|) that has the highest critical temperature be- 
low Xnmp, see Fig. ^ One can show that as r varies from 
+00 (no dipole interaction, Tt — 0) to (strong dipole 
interaction, X t = X s ), Xnmp moves from Xloff towards 
T s . The parameter p is the measure of the triplet compo- 
nent admixture, which is zero at X = Xnmp and increases 
as temperature decreases [for the parameter values used 
above, p(T = 0) ~ 0.23]. 

In order to determine the spatial structure of the NMP 
phase just below the critical temperature, one needs to 
evaluate the fourth-order terms in the free energy density: 

X 4 =/3 1 (|r 7l | 4 + |r; 2 | 4 )+/3 2 |r, 1 | 2 |r ?2 | 2 . (12) 

The coefficients /3i 2 are functions of temperature. The 
above expression is positive definite if (3\ > 0, /3 2 > — 2/?i. 
At 02 > the minimum is achieved in the state (771 , ?y 2 ) ~ 
(1,0) or (0, 1), while at (i 2 < one has (771,772) ~ (1, e^), 
where ip is an arbitrary phase. Near Xnmp one can set 
p = q c = in the expressions for f3\ 2 , which gives 0i = 
/? 2 /4 = -(AM^)/32^ 2 X2 Mp )Re#' ; (l/2 - tz NMP ) > 0. 
Therefore, the NMP phase transition in the vicinity of 
Xnmp is of the second order, with (ip,d z ) ~ (l,p)e lqcZ . 
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The determination of the full phase diagram, including 
the spatial structure of the order parameter at all T < 
Tnmp, is much more difficult, because of the competition 
between the NMP, LOFF, and uniform BCS phases. We 
leave this problem for future studies. 

The origin of the NMP instability can be elucidated 
using the Ginzburg-Landau expansion of the free energy 
|jHJ at small h and q = qz in the vicinity of T s . In addi- 
tion to the usual uniform and gradient terms, the energy 
density contains mixed-parity terms, which are linear in 
gradients 



ip*A s ip + d* z A t d z + iKh(ip*V z d z + dy z ip), (13) 



where A a = N F {T - T a )/T a - K a V 2 z , 
K a = 7C(3)iV F (^ u 2 )/167r 2 T 2 ) and K = 

7((3)N F ((j) s <j)tV z } /4tt 2 T s 2 . Because of the last term, 
which favors a spatial modulation of the order pa- 
rameter, at h > h NMP = yj NpK s (T s - T t )/T t K 2 the 
superfluid phase transition is of the NMP type. In con- 
trast, the onset of the singlet LOFF instability is marked 
by the coefficient K s changing sign at /iloff > /inmp- 

To summarize, we found that in a homogeneous Fermi 
gas with s-wave attraction, the atomic dipole interaction 
produces triplet pairs that can mix with the singlet order 
parameter in a spatially modulated superfluid state. If 
the densities of atoms in different hyperfine states are un- 
equal then the nonuniform mixed-parity state has higher 
critical temperature than the LOFF state, even at weak 
dipole interaction. To relate these findings with the cold 
atom experiments the effects of the trapping potential 
should be included. Another interesting open question 
concerns the fate of the NMP state in the strong-coupling 
regime near the BCS-BEC crossover. 

Finally, we would like to remark that the possibility 
of a nonuniform singlet-triplet mixing is not restricted to 
polarized Fermi gases. For example, in a fully isotropic 
system the pairing interaction V(k, k') in Eq. (J2J de- 
pends only on the angle between k and k' and can there- 
fore be factorized in terms of the spherical harmonics 
Y; m (k). In the singlet s-wave channel the order param- 
eter is a scalar tjj, with </> s (k) = Ybo(k) = 1, while in 
the triplet p-wave channel the order parameter is now 
represented by three vectors d m , with m = 0, ±1, corre- 
sponding to (f>t t m(k) = Yi m (k). When expressed in terms 
of the Cartesian components of k, the triplet order pa- 
rameter becomes a 3 x 3 complex matrix Ai a |20| . In the 
presence of magnetic field h, there is a contribution to 



the free energy density of the form ih^ifj*^ a Ai a — c.c), 
which is invariant under all required symmetry opera- 
tions: orbital and spin rotations, time reversal, and in- 
version. Similarly to the last term in Eq. (|13(l this can 
lead to the NMP instability at sufficiently high h. 
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